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MnAs/GaAs superlattices, made by 5-doping GaAs with Mn, are known as digital ferromagnetic 
heterostructures. Here we present a theoretical density functional study of the electronic, magnetic 
and transport properties of such heterostructures. In the absence of intrinsic donors these systems 
show an half metallic density of states, with an exchange interaction much stronger than that of a 
random alloy with the same Mn concentration. Ab initio ballistic transport calculations show that 
the carriers with energies close to the Fermi energy are strongly confined within a few monolayers 
around the MnAs plane. This strong confinement is responsible for the large exchange coupling. 
Therefore the system can be described as a two dimensional half metal with large conductance in 
the MnAs plane and small conductance in the perpendicular direction. 

PACS numbers: 75.50.Pp, 71. 15. Mb, 73.23.Ad, 73.50.-h 



Diluted magnetic semiconductors (DMS) based on III- 
V materials |J § are receiving great attention for the in- 
triguing interplay between magnetic and electronic prop- 
erties H and for the large potential for future electronic 
devices W. Gai-xMn^As, which is structurally compat- 
ible with most epitaxially grown III-V's, is the prototype 
of these materials. Here the Mn ions substitute randomly 
for the Ga ions, providing both localized 5 = 5/2 spins 
and free carriers (holes) . The antiferromagnetic exchange 
coupling between the hole spin and the Mn spin is respon- 
sible for the long range ferromagnetic interaction, which 
can be described in a first approximation by the Zener 
model H . One important result of the model is that the 
Curie temperature is a linear function of the Mn concen- 
tration. Although the inclusion in the model of correla- 
tion effects Q modifies this linear dependence, there is 
general agreement that an increase in the Mn concentra- 
tion enhances the Curie temperature of the system. 

Unfortunately the solubility limit of Mn in GaAs is 
rather small, since MnAs naturally occurs in a hexago- 
nal phase with [6]-coordinated Mn [Q. However large Mn 
concentration, up to 50% can be obtained in zincblende 
MnAs sub-monolayers embedded into GaAs to form a 
MnAs/GaAs super lattice ||. These structures, called 
digital ferromagnetic heterostructures (DFH) have re- 
markable properties. First, the Curie temperature decays 
for increasing GaAs inter-layer thickness, and saturates 
for thicknesses larger than ~50 GaAs monolayers. The 
value of T c at saturation depends on the Mn concentra- 
tion A saturation is unexpected according to the 
Zener model for three dimensional systems, since the to- 
tal Mn concentration in the sample decreases with the 
increase of the GaAs thickness. This separation depen- 
dence suggests that DFH's behave like planar systems, 
although to date there is no direct proof that the free car- 
riers responsible for the magnetic interaction are confined 
in the MnAs layers. Secondly, Hall measurements in the 



direction parallel to the MnAs planes show an anomalous 
Hall effect for undoped samples, which disappears upon 
Be-doping j|, |l0). Large Shubnikov de Haas oscillations 
are found in doped samples, although surprisingly the 
charge densities extracted from the Hall coefficient and 
from the Shubnikov de Haas oscillations are different. 
This suggests that two different carrier types could be 
present in the system. 

In this letter we investigate theoretically the magnetic 
and transport properties of undoped DFH. In particu- 
lar we address the following questions: i) what is the 
real dimensionality of the system? ii) are the carriers 
spin-polarized? iii) what is the carrier distribution in the 
system? 

We perform density functional theory Jll[] calculations 
in the local spin density approximation. The numeri- 
cal implementation included in the code siesta 12 uses 
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pseudopotentials and a localized pseudoatomic orbital 
basis set, and allows us to investigate diluted systems 
within a supercell approach. Extensive details of the 
method and its optimization have already been presented 
p3| . Each DFH superlattice is constructed by consid- 
ering N GaAs cubic cells (8 atoms in the cell) aligned 
along the z direction (see figure |]a). One Ga plane is 
substituted with Mn and periodic boundary conditions 
are applied. This leads to an infinite MnAsi/GaAs27v-i 
superlattice, where MnAs zincblende monolayers are sep- 
arated by a 5.65 x NA thick GaAs layer. 

In figure [l] we present the band structure for the case 
N = 8 (with a 45. 2A thick GaAs inter-layer), for both 
the majority and minority spins. The points X\ and 
X2 are in the MnAs plane and denote directions along 
the edge and the diagonal of the cubic supercell respec- 
tively, r — > Z\ is the direction orthogonal to the MnAs 
plane. The band structure of figure |lj is that of an half 
metal, with a gap of approximately 0.65 eV in the minor- 
ity spin band. The majority spin band presents a strongly 
anisotropic dispersion. It is metallic and quite dispersive 
in the MnAs plane (directions X\ — > V — > X2), while in 
the direction orthogonal to the MnAs plane there is the 
formation of a narrow impurity band. Since the Fermi 
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FIG. 1: Band structure for a MnAsi/GaAsis superlattice: (a) 
majority and (b) minority spins. X\ and X2 are in the MnAs 
plane and denote directions along respectively the edge and 
the diagonal of the cubic supercell. V — > Z\ is the direction 
orthogonal to the MnAs plane. The horizontal line indicates 
the Fermi energy. 



energy Ep cuts through the impurity band, MnAs/GaAs 
looks like a two dimensional half metal with small hop- 
ping between the MnAs planes. The width of the impu- 
rity band decreases with increasing separation between 
the MnAs planes (665 meV, 412 meV and 256 meV re- 
spectively for N=4, N=6 and N=8), which is consistent 
with the decrease of the overlap between the wavefunc- 
tions centered on the MnAs planes jl4| . 

We calculate the strength of the magnetic coupling by 
calculating the difference Afa between the total energy 
of supercells with either antiferromagnetic or ferromag- 
netic alignment of the Mn ions [ p| . We find that the 
ferromagnetic interaction is much stronger than that in 
a random alloy with the same Mn concentration. For 
instance in the MnAsi/GaAsi5 DFH with a Mn con- 
centration of cc=0.06, Afa is 515 meV, whereas in the 
corresponding random alloy it is only 160 meV. More- 
over, Afa seems to be rather independent of the GaAs 
thickness for the range of thicknesses investigated here 
(531 meV, 533 meV and 515 meV respectively for N=A, 
N=6 and A^=8), suggesting that the carriers responsible 
for the ferromagnetism are indeed strongly confined in 
the MnAs planes. 

We also investigate the effect of introducing As anti- 
sites (AsQa) by considering a MnAsi/GaAsn superlat- 
tice with half a monolayer of Aso a close to the MnAs 
plane. This corresponds to one Asc a for two Mn, hence 
according to the nominal valence, to compensation. In 
contrast to what expected from the Zener model, the 
coupling between the Mn ions is still ferromagnetic with 
Afa=70 meV. This result is similar to that found for 
the tetrahedral GaiMn2As2 complex in GaAs |l5| , al- 
though in that case Afa=20 meV. We interpret the fer- 
romagnetism at compensation in presence of Asc a as 



due to localized carriers Zener-coupled to the Mn ions 
|l5| . The MnAs/GaAs superlattices have stronger ferro- 
magnetic coupling than that of the GaiMii2As2 complex 
since the carrier localization close to the MnAs plane is 
very strong. These results suggest that also in DFH the 
Curie temperature is determined by a delicate interplay 
between Mn ions and intrinsic defects. 

We now move to the transport properties. Here we 
generalize the technique of reference [[l6| to the case of 
non-orthogonal tight-binding model with singular cou- 
pling matrices. The conductance is calculated in the bal- 
listic limit with both the current in the Mn plane (CIP) 
and perpendicular to the Mn plane (CPP). We first ex- 
tract the tight-binding Hamiltonian and overlap matrix 
giving the correct charge density. The matrix elements 
are computed by numerical integration over a real space 
grid [ p^j . The resulting tight-binding Hamiltonian com- 
prises in principle an arbitrary large number of nearest 
neighbor interactions, which is fixed by the cutoff ra- 
dius of the pseudoatomic orbitals. We then rewrite both 
the Hamiltonian and the overlap matrix in a tridiagonal 
form. This procedure fixes a "natural" length scale. In 
fact we redefine the unit cell in such a way that each unit 
cell is coupled only to nearest neighboring cells. Hence 
the problem is reduced to an effective quasi-one dimen- 
sional problem. The scattering matrix is then calculated 
by a well established Green's function technique Jl6| , and 
the spin conductance by using the Landauer-Biittiker 
formula jlTj and integrating over the two dimensional 
Brillouin zone in the plane orthogonal to the direction of 
transport, 



r CT = y]TTri ff (fc||)^(fc||)t 



(1) 



Here fc|| spans the two dimensional Brillouin zone in the 
plane orthogonal to the direction of transport, and t a (k\\) 
is the fc|j -dependent transmission matrix for the spin a. 
We assume a two spin fluid model where there is no mix- 
ing between majority and minority spins. It is important 
to point out that in the CPP case our supercell geometry 
leads to narrow bands. It is then crucial to perform large 
/c|| -point sampling. Here we consider up to 800 fcii -points 
in the two dimensional irreducible Brillouin zone. 

In figure [| we present the conductance per unit area 
as a function of the position of the Fermi energy for a 
MnAsi/GaAsi 5 superlattice for both the CIP and CPP 
directions and both spins. We also project the conduc- 
tance onto the atomic orbital basis set in order to deter- 
mine the orbital character of the electrons carrying the 
current |H. Moreover for the CIP configuration in fig- 
ure we show the corresponding spin polarization of the 



conductance £, which is defined as £ = 



The cou- 



rt +n 

ductance shows an half metallic behavior for both the 
CIP and CPP directions. This means that there is a 
100% spin polarization of the conductance for undoped 
samples (Ep is in the gap for the minority spin). In fig- 
ure p| we also show the position of the GaAs valence band 
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FIG. 2: Total and partial conductance per unit area for a 
MnAsi /GaAsis superlattice as a function of the position of 
the Fermi energy for majority ( j) and minority (J.) spin bands. 
The vertical line denotes the position of the Fermi energy for 
undoped samples. 
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FIG. 3: Spin polarization of the conductance for the CIP 
configuration. The vertical lines denote the position of the 
Fermi energy (solid line) and of the GaAs valence band edge 
(dotted line). 



edge, which is obtained by aligning the lower lying As s 
levels of bulk GaAs and the MnAsi/GaAsis superlattice. 
It is easy to note that £ is 100% for energies close to the 
GaAs valence band edge and that the polarization per- 
sists well within the valence band. For this reason DFH 
appear as good candidates for spin injection into p-doped 
GaAs. 

If we now analyze closely the conductance of the ma- 
jority spin band we notice that it is strongly anisotropic, 
as expected from the bandstructure of figure [l]. In the 
CIP direction the conductance as a function of the po- 
sition of the Fermi energy is rather flat with an orbital 
contribution mainly from the As p orbital. There is also 
a significant contribution (roughly 20% of the total con- 
ductance) coming from the Mn d orbitals, suggesting 
that the carriers in undoped samples are confined close to 
the MnAs planes. It is also important to point out that in 



fc-space the contributions to the conductance are rather 
uniformly distributed over the two dimensional Brillouin 
zone in the plane orthogonal to the transport direction. 
These are strong indications that MnAs/GaAs is a good 
isotropic metal in the MnAs planes. 

For the CPP direction the situation is rather differ- 
ent. The conductance as a function of the position of Ep 
shows clear one dimensional features and for undoped 
samples it is solely given by charge in the impurity band. 
The orbital content of the impurity band is mainly As p 
with small contributions from Mn d t2- The conductance 
at Ep comes from few fc-points around the T point. In the 
case of carriers with parabolic dispersion this happens for 
scattering at a potential barrier, when the carrier energy 
is very close to the top of the barrier. In fact carriers 
with largest kinetic energy in the direction orthogonal to 
the barrier (small k in the parallel direction) have largest 
transmission amplitude. This situation is similar to that 
of tunneling junctions and so we describe the transport 
as tunneling-like. 

In summary, from ballistic transport calculations, 
MnAs/GaAs appears to be a two dimensional half metal 
in the MnAs plane, with tunneling-like conductance be- 
tween the MnAs planes. In order to better understand 
the spatial arrangement of the current in figure ^| we 
present the charge density distribution in real space, 
p(r), calculated only for those states contributing to 
the conductance and within a 0.6 eV range around Ep 
{E F - 0.3eV < E < E F + 0.3eV). We also notice that 
the individual contributions from occupied (Ep — 0.3eV < 
E < E F ) and unoccupied (Ep < E < E F + 0.3eV) states 
are very similar. 




FIG. 4: Charge density distribution in real space for 
MnAsi /GaAsis (a) calculated only for those states contribut- 
ing to the conductance and energy within 0.3 eV from Ef: 
CIP (b) and CPP (c) configurations. 

The figure confirms that the current in the CIP case is 
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distributed mainly in a narrow region around the MnAs 
planes, with small spillage outside. Planar averaging Q 
of the charge density of figure || shows that the charge 
density reaches ~0.1 of its MnAs plane value within only 
three GaAs monolayers from the MnAs plane. In con- 
trast, the CPP current is mainly located at the Mn plane 
with small contributions from the GaAs layers. This 
means that carriers are strongly confined in the MnAs 
plane and the perpendicular transport is via hopping be- 
tween the planes. 

Finally in figure ^| we present the macroscopic average 
]18[ of the Hartree potential. The macroscopic average 
is obtained by first taking a planar average and then by 
averaging the result over the GaAs lattice spacing. This 
does not correspond exactly to the total potential which 
the quasiparticles feel in the structure, but it gives a good 
indication of the latter |H . The potential profile shows a 
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FIG. 5: Macroscopic average of the Hartree potential for 
MnAsi/GaAsis. 

deep well in the MnAs plane separated from a flat region 
by two barriers located one monolayer from the well. The 
presence of the barrier suggests that n-type doping in the 
GaAs region will only weakly affect the MnAs planes, 
since the extra electrons need to tunnel through a large 
barrier to access the MnAs planes. 

In conclusion we have studied with ab initio methods 
the magnetic and transport properties of GaAs/MnAs 
digital ferromagnetic heterostructures. We have shown 
that they are two dimensional half metallic systems with 
large metallic conductance in the MnAs plane, and small 
tunncling-like conductance perpendicular to the MnAs 
plane. We have shown that the electronic charge at the 
Fermi energy is mainly distributed within a few monolay- 
ers around the MnAs plane and that the Hartree poten- 
tial shows a double barrier profile. This potential profile 
can sustain the presence of a two carrier population for 
DFH's that are Be-doped in the GaAs region. 
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